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Abstract 

Using a contour integral representation we analyze the multi-instanton sector 
in two dimensional U{N) Yang-Mills theory on a sphere and argue the role of 
multi-instanton in the large phase transition. In the strong coupling region 
at the large A'^ , we encounter "singular saddle point". Because of this situa- 
tion, "neutral" configurations of the multi-instanton are dominant in this region. 
Based on the "neutral" multi-instanton approximation we numerically calculate 
the multi-instanton amplitude , the free energies and the Wilson loops for finite 

. We also compare our results with the large exact solution of the free 
energy and the Wilson loop and argue some problems. We find the "neutral" 
multi-instanton contribution bridges the gap between weak and strong coupling 
phase. 



1 Introduction 



In the past years, various attempts to solve low energy QCD have done. Among them 
key idea is a kind of string description of QCD because the string description automat- 
ically implies the confinement. Toward this end , there was remarkable progress for 
QCD in two dimension in the last few years. In particular Gross and Taylor discovered 
the string descripton of two dimensional SU{N) or U{N) Yang-Mills theory in terms of 
expansion It is well known that the partiotion function ( and the Wilson loop) 
of the Yang- Mills theory on genus G Riemann surface Eg can be expressed by the form 
of the sum over irreducible representations of the gauge group 0. The form is called 
the heat kernel representation . They showed that expanding the heat kernel by , 
one gets the sum over the branched covering maps of the Riemann surface Eg weighted 
by A^^~^^ exp(— lAcover)- Here g and Acovev are genus and area of the covering space 
respectively. That is, a no- fold string theory with target space Sg. But later it was 
discover by Douglas and Kazakov that on a sphere this system undergoes a 3rd order 
phase transition in the large N limit 0. They showed whereas in the strong coupling 
phase the large N solution of the free energy has the string expansion corresponded 
to the map from sphere to sphere , in the weak coupling phase the large solution is 
trivial and has no such expansion. Thus it is impossible to interpret QCD as the string 
theory in the weak coupling phase. 

Viewed from the strong coupling phase ,the string sum becomes divergent at the 
phase transition point due to the entropy of the branch points^. So the entropy of the 
branch points physically causes the phase transition. On the other hand , viewed from 
the weak coupling phase instanton induces the phase transition|^, |]. We recall Witten's 
work on two dimensional QCD0. He showed that the Yang-Mills partition function on 
Riemann surface can be expressed by a sum over the instanton. The weight associated 
with the instanton is corectly determined from the heat kernel representation!^. Using 
this formula Gross-Matytsin showed the following results[§]. The instanton sector 
gives the weak coupling result fot the free energy. In the weak coupling region the 1 
instanton amplitude of charge ±1 is strongly suppresed in the large limit but at the 
phase transition point the damping factor disappears . 

In addition one can see using their results whereas the damping factor is disap- 
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pearded at the phase transition point, the 1 instanton amphtude of charge ±1 still has 
order A^~^/^ in the strong coupling region. And it's contribution to the free energy is 
order N~^^'^ and thus negligible in the large N limit. Hence the 1-instanton effect is in- 
sufficient to explain the large N phase transition and it is important to investigate the 
effect of the multi-instanton in the large N limit. That is our motivation for analyzing 
the multi-instanton sector. 

Since in the strong coupling phase Yang -Mills theory has the string description, by 
studying the phase transition we can see how the nonperturbative effect of the multi- 
instanton construct the string picture. That is instructicve for real four dimensional 
QCD. In four dimensional QCD it is unlikely that there is such phase transition|0. 
But perturbation can not lead to the string picture, it is important to know how 
the nonperturbative effect construct the string picture. By considering the case of 
two dimension in detail, there might give some insight for four dimensional QCD. We 
should keep the above thing in mind when we study the large N phase transition in 
two dimensional QCD. 

In two dimension the Yang-Mills theory has no transverse gluon and it is in some 
sense topological. So coupling to various matter fields is very important. But such sys- 
tems are obtained by considering random walk for the Wilson loop of the Yang-Mills 
theory]^, our method might become another approach to the above systems. 



In the preceeding paper the auther gave a scenario of the large phase transition 
which explains the 0(1) contribution to the free energy in the strong coupling region 
using the multi- instanton 0. In this paper following the preceeding paper, we continue 
our analysis on the multi-instanton contribution to the free enegy and the Wilson loop 
and find out that there is a sort of neutrality in the strong coupling phase at the 
large A^. That is, "neutral" configurations of the multi-instanton are dominated for 
both the free energy and the Wilson loop in this region. Based on the "neutral" multi- 
instanton scheme, we also calculate the free energy and the Wilson loop for gauge group 
[/(3), [/(4)and U{5). And we compare our results with the large A^ exact solution of 
the free energy and the Wilson loop|TD[ and find that the neutarl multi-instanton 
bridges the gap between the weak and the strong coupling phase in two dimensional 
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QCD on sphere. 

The content of this paper is as follows. In section 2 the formula of the multi- 
instanton amplitude and the partition function are obtained in terms of contour inte- 
gral. In section 3 the formula of the Wilson loop average are also obtained in terms 
of contour integral. In section 4 nature of the 1-instanton and the dilute gas approx- 
imation are considered. In section 5 the structure of the multi-instanton amplitude 
in the large limit is analyzed in detail. In section 6 numerical calculation of the 
multi-instanton amplitude , the free energy and the Wilson loop are performed. Also 
our results are comared to the large N exact solution . In conclusion we sumary the 
results. 



2 partition function 

The partition function of two dimensional U{N) Yang-Mills theory has so called the 
heat kernel representation. The representation was first considered by A.Migdal on 
disk in the context of the real space renormalization group and later generalized to 
arbitraly 2-manifold by B.Rusakov0. On a sphere with area A , the partition function 
becomes 

ZiA) = JvA^expi-^Jd'x^tTF,,Fn 

= E(dim/2)^exp(-^C2(i?)). (1) 

R 

Here R is irreducible representation of the gauge group U{N),dimR is the dimension 
of representation R and C2{R) is the value of the second Casimir operator of rep R. In 
terms of the highest weight components of R , rii > n2 >■■■> nj\f , dimi? and C2{R) 
are given by 

dimi? = n (1-^^), (2) 

l<i<j<N ^ J 

N 

C2{R) = 5]n,(ni + Ar + l-22). (3) 
1=1 
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If we define k as Ui — i -\ — ^ , the partition function becomes [| 

Z(A) = const e^(^^-^) ^ U A2(/) exp(--^ ^ ^D- (4) 

We remark that the condition that /j's satisfy the relation /i > /2 > . . . > /at is irrevant 
because the configuration such as /j = Ij does not contribute to the Van der Monde 
determinant A and Z{A) has the index permutation symmetry. Then the sum becomes 
free sum over for odd or sum over (Z + for even N and we can rewrite it by 
the Poisson resumation formula. This is a sort of duahty transformation. We get the 
following expression f\ 0: 

Z{A) = e^(^^-i)(^)^^ Y: eSti-e-^St.™?^({m}), (5) 

N 



n dy^ 11(4 - 47r2mJ)e-^ ^» 

= / Y[dy^A{y, + 27^m,)A{yi-2nmi)e-^^^=ly' (6) 
•^-°°i=i 

where yij = yi — yj, rriij = rrii — rrij, e = 1 for odd and e = — 1 for even N. In 
P, H, 01 it was showed that {m} corresponde to all Euclidean classical solutions (which 
we call instanton) up to gauge transformation . So the instanton is dual to the highest 
weight component. The instanton have nonperturbative effect with respect to both 
and A. Hereafter we call instanton charge and call number of nonzero mj's 
instanton number. In the following we analyze the phase transition using the duality 
and see the instanton induced the phase transition. That is reminiscent of the duality 
transformation of classical XY model and the Kosterliz-Thouless phase transition . 
In this case viewed from the low temperature the vortex condensation causes the phase 
transition. 

The multi- instanton amplitude w({m}) looks like the partition function of the Gaus- 
sian Hermite matrix model but there is a deformation in the Van der Monde determi- 
nant A. In this section we rewrite w{{m}) using the method of ortho-polynomial (in 

^Hereafter we absorb A into A . 

similar expression was obtained by M. Caselle et al jl^. They showed that the phase transition 
is due to the winding modes of "fermion on circle" . 
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this case Hermite polynomial) and obtain a new contour integral representation. The 
new representation makes clear the role of the multi-instanton in the large N phase 
transition. 

First using the property of the Van der Monde determinant we obtain : 
w{{m})= sgn^ H dyiPa(i){yi + 27rmi)Pi^„a(i){yi-27:mi)e 2^^S (7) 



where Sn is the permutation group on N object and 



1 rr . N_ 



= r,r,r NM ^niy TTJ^) = ^" + (lower power) (8) 



is the ortho-polynomial under the Gaussian weight ; 

J —oo i V 

In the following let us assume number of non-zero m^'s is k i.e. number of instanton 
is k. Using the permutation symmetry nonzero m^'s are driven to m^'s from i = 1 to 
i — k. Then we get, 

w(mi,...,mfc,0, ...,0) 

^ ^ N 2 

= ^ sgn/i ^ Yl ^cT(i)(^a(i),Moa(i) n / dyiP^^i){yi + 2TTmi)Pi^^^(i){yi-2TTmi)e'^yi 
= {N-k)\J2^SW ho- ■ -hai- ■ -ha,^- ■ -hN-i 

^ f°° N 2 

X H / dyiPaXVi + 27rmi)P^(„.)(yi - 27rmj)e-5Afi . (10) 

Here /i & Sk is the element of the permutation group acting on the set {oi, . . . , a^}. 
Using the Taylor series expansion of the Hermite polynomial, we obtain for the above 
integral, 

^ /"OO ^ 2 

n / dyiPa^iyi + 27rmi)P^(„.)(yi - 27rmi)e'^yi x • • • /i„J-^ 
= TT V /^"'^ ^' (^^^)°'-^'(-27rm)^(°')-\ (11) 
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The each scries in the above equation can be represented as a contour integral by the 
following transformation formula. 

min(a,6) ^ r rif 1 ^ 

E 7:^i^rw-' - <f &,^"'7(irm + (12) 

i=0 



where the contour of t encircles only the origin counterclockwise. Using this formula, 
we obtain the following contour integral representation for the /c-instanton amplitude. 

w(mi, . . . ,mfc,0, . . . ,0) 

{N-k)\^ .N. ^1 _ rdh fdtk 1 



mi ti ruk tk 

Here ZciP) = / H dxA'^{x) exp(— | J2 x'^) is the partition function of the Gaussian Her- 
mite matrix model. We remark the configurations such as ai = 02 do not affect the 
above equation. Hence we can replace Z]^^^..^afc independent sum I]^^^ ■ ■ -S]^"^- 
In that expression there are many pole free terms. After eliminating them, the ampli- 
tude becomes 



iN-k)\ N rdt^__ r ^,-^(mit.+-+n.it,)n + 1)^ . . . (i + 1) 



N 



X X] sgn/x ^ ■ ■ ■ m . (14) 

The last sum is just determinant of matrix M which has (ij) element 

= — 3— . (15) 

It's diagonal element is 1 and this determinant has all information about the interaction 
between instantons. 

For the sake of numerical calculation we further rewrite the /c-instanton amplitude 
in terms of another contour integral. For this purpose, we classify elements of the 
permutation group 5"^ into the conjugacy classes. By symmetry argument, the contri- 
bution from different elements belonging to same conjugacy class is same. Then we 
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obtain, 



J2 e"''+-+""^e-^^"''+-+<^w{mi, . . . , m^, 0, . . . , 0) 

mi,...,mfc^O 



N\ A' ^ ,^ J 2m J 2m 

X sgn[(7]r[(7]Mi^(i)M2<,(2)---Mfe^(fe), (16) 

conj. class 

where 

= + — ^ ^ — ^(t + -) - log(l + -). (17) 

We assume a has cycle structure [l"'^ • • • A;'^*=]((Ti + • • • + kuk — k) and put T[a] as 
number of elements in the conjugacy class which a belongs to, 

sgna = (-)'^^+^^+-+%$i, (19) 
E = E • (20) 

conjclass o-i , ■ ■ ■ .^^fc^O 

cr^H |-fccr^ = fc 

We use (12345) (67) type element as the representative element of the conjugacy class. 
Because the multiple integral factorize according to the cycle structure of cr, the multi- 
instanton amplitude is exponentiated with a constraint (cti + • • • + kuk — k) which is 
expressed by a contour integral of z. We get the following result: 

^ . . . , mfe, 0, . . . , 0) 

mi,...,mfe^O 

= iNCk^Zai^) / cxp(f: ^^z^a,{A)), (21) 

where the contour of z encircles only the origin counterclockwise and 

Tni,...,mj^O 

(22) 

is the "connected" amplitude of the j instantons. In the eq(21) we extrapolate the 
range of summation inside the exponential from J2'j=i to J2jLi jbut this does not affect 
the equation due to the z integral. 
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Using the above equation we finally obtain for the partition function 

k=0 mi ,...,mfc7^0 

= ^wcak(A)/— p^^exp(E^^^«,(A)), (23) 

j=i J 

where 

^weak(^) = e^(^'-i)ZG(A/iV) = conste^'(^-^^°s^)-^. (24) 

Here we use the the dual property of the Hermite Gaussian matrix model as Zg{(3) = 
p-N^ Zq{1/ j3). Since Z„cak(^) is the 0-instanton (i.e. all mj = )contribution to the 
partition function , Z„eak(^) is just the equation (4) before the duality transformation 
in which the discrete sum is replaced by the continious integral. That is the peculiar 
property of the Poisson resumation formula. Douglas-Kazaov showed that Z„eak(^) 
is in fact the partition function in the weak coupling phase at the large A^. Thus the 
multi-instanton effect is essential in the strong coupling region at the large N. 



3 Wilson loop 

In this section, following the preceeding method , we give the contour integral represen- 
tation of the multi-instanton sector in the Wilson loop. The heat kernel represenatation 
gives the following expression for the Wilson loop on a sphere : 

= ^Edimi^dim^e- "^'^^^"'^^^"^'^^''' ^ / dUxR{U)xs{U^)tTU\ (25) 

Here the loop devide the sphere with area A into two disks with area Ai,A2. The 
above holonomy integral gives the linear combination of the Clebsh-Gordan coefficient. 
For example of n = 1 case , The sum over representation R, S becomes that over R 
multiplied by the sum over all posible attachment of one box to the Young tableau of 
R. Hence in the same way as the partition function , the sum becomes free sum over 
for odd iV or (Z + i)^ for even N. By the Poisson duality we get the following 
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result of the Wilson loop average 
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Wn{Ar,A2)Z{A) 



6 24 



A N 



X 



J — CO iV iV 



where Hk=i exp(— 27rmmfc^) is the classical value of the Wilson loop with the multi- 
instanton configuration. From the above equation, the Wilson loop average apparently 
has Ai ^ A2 exchange symmetry. We also find that calculating the Wilson loop is like 
calculating the partition function with additional imaginary charged instanton. Hence 
in the same way as the partition function, we can rewrite this in terms of a contour 
integral. In the /-instanton sector ( i.e. number of non-zero m's is / ) , we devide the 
place where the imaginary chage put on . One is k with nonzero rrifc, another is k with 
zero mfc. Hence the above equation becomes 

Ao r°° IL. N v^iV 2 

x[/e-2™^^ / []%e-5AE,=i2^« 

X A(yi + 27rmi + mAg/A^, y2 + 27rm2, ■■■ ,yi + 27imi, yi+i, ■ ■ ■ , yjv) 
X A(?/i - 27rmi + inAi/N, y2 - 27rm2, ■■■ ,yi- 27imi, yi+i, ■ ■ ■ , i/tv) 

+ {N-l) / J^dy^e-^S-i^' 

X A(yi + 27rmi, ■■■,yi + 2™^, yi+i + m^a/A^, yi+2, ■■■,yN) 

xA{yi - 2'Kmi, ■■■,yi- 27imi,yi+i + inAi/N,yi+2, ■ ■ ■,1/Af)]- (27) 

The first integral in the big parentheses, using the preceeding method, has the following 
contour integral representation : 

' ^^-Zai-) / ^ ■ ■ • / ^e-^K*.+-+"^^0(i + 1)^ . . . (1 + 



N\ ' A' J 2ni J 21X1 ti 

xe2™i^^*i-^^*Met(')5. (2^ 
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The second integral also has the contour integral representation as 

~ ^ ~ Zg(-) / ^ • • • / ^e-^K*^+-W*0(i + l)iv . . . (1 + 1) 
A^! J 27ri / 2ni ti 



N 



where 



xe-^^*'+^(l + J-)^det('+i)C, (29) 

Bii — T , bi — — 27rmj + ■^^rr-'^ii, (30) 

|^(l + t,)-ti ^ TV ^ ^ 

^^J' ^ cj^-^^"^^ _^ , Cj = -27rmi(l - 6i,i+i) + ^-^^S^+i. (31) 

For the sake of numerical calculation , we further rewrite it to a exponential form. The 
first term of the /-instanton contribution to the Wilson loop has the following form: 

J2 ^rm+-+m, I ^e-iv5^,(*0/(mi, ii)det«S. (32) 

mi,-,mi^Q '' i=l 

In the determinant we assume the cycle including the index 1 to have length k and 
explicitly to be (1, ai, . . . , ak-i) .Thus the determinant will be 

det(^)S = ^ • • • 5a,_iidet('-')5, (33) 

fe=l ai^---afc_i^l 

here det^'"''^^ is the determinant not including the index 1, oi, . . . afc_i. Under the 
index rename invariant measure , we can rename Oi, . . . , au-i to be 2, . . . , /c. Thus wc 
get 

det«5 ^ 7!^^(-l)'"'^i2i^23 • • • BkidS'-^^B (34) 
Hence eq.(32) becomes 

±r dt 



•••,m,^0 i=l 

/!E^A-fe(^,^) E e-+-+-^/n?:e-^*'^«(*')det(^)M. (35) 



k=0 ' mi,---,mk^O 

where 

^ dt, 



--, . Ao , „ . Ai —An , 7,2 Ai Ao , 



10 



is the type 1 connected amplitude of /c-instanton 

In the same way , under the index rename invariant measure ,we get 
i+i /I 

det('+i)C ^^^(-i)'=-iCi2C23- ■ •G+ndet('+i-'=)a (37) 

k=i [i-k + iy. 

Hence the second term of the Z-instanton contribution to the Wilson loop becomes 

J2 / Yl ^^-N^^AU) I ^e-'^*'+^(l + ^)^det('+^)C 

,„ J A 1 27ri J 2,TTi ti+i 

= ^!Er[7z-.(^i,^2) E e-^+-+-^/n|^e-^*™'(*»)det('^)M. (38) 

k=0 ' mi,---,mf.j^O i=l 

where 

1 ^ 7TZ -I 

mi,---,mfc^u 2=1 «=i 

xe-=^*'=+^ (1 + J_fc^2C23 ■ ■ ■ Cfc+n. (39) 

is the type 2 connected amplitude of k-instanton. Using the fact that /c-instanton 
amplitude can be rewritten as eq (21), we finaly obtain the following result: 

WniA,, A2)Z(A)/Z^eak(A)e-^^ 

= ir^Ef^(L M + E 3TI )exp(E— — ^'«i(^))(40) 

This is the base of our numerical calculation of the Wilson loop average. 

For example , in the 0-instanton sector only 70 contributes to the Wilson loop. We 

get 

W^„(Ai, A2)Z(A)/Z„,,k(A)e-^|o-inst = ^70 = ^ / Ae-^*(1 + V. (41) 

i\ IS J 2m t 

By rescaling t — > Nt, we obtain the known result of the Wilson loop in the weak 
coupling phase at the large N [l^ : 



W^.(Ai,A2)|o_inst = ^^^Ji(2n^^) + 0(l/iV). (42) 

So the multi-instanton effect is essential for the Wilson loop in the strong coupling 
phase. Apparently if we extrapolate the above equation to A — 00 ^ with fixed Ai , 
we cannot obtain the simple area law on infinite disk but obtain a oscilating behavior. 

■^Thc phase transition occures before A ^ 00. 
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4 Effects of 1-instanton 



In this section we consider the contribution of one instanton to the partition function]^ 
and the Wilson loop. We see that the 1-instanton effect is insufficient to explain the 
large N phase transition. We use the word "1-instanton" to mean that only ai, f3i, 70, 71 
are included in the partition function and the Wilson loop. So this is a kind of dilute 
gas approximation in our system. 

The dilute gas approximation leads to the following result for the partition function. 

CT) ^ f 2^i^r^^^p^^"^(^)) - ^''^ 

We should note since maximum number of the instanton is , the partition function 
in this approximation is not completely exponentiated as exp(Q;i). In the large limit, 
we can apply the saddle point method to ai. The saddle point equation = has 

two solutions 

-1 ± a/i ^ 



t = t±^ . (44) 

We see quite different behavior according to A^vr^m^. 

For the case of A < vr^m^, both saddle points exist on negative real axis and satisfy 

$„(t_) < ^Ut+). (45) 

Neverthless only saddle point is selected since t- steepest decent line is not consistent 
with the original contour. Hence we obtain in the large N limit , 



^ '••/27riV|<l>^^^(t+)| 



(46) 



where 7(x) is defined as 



X 1 + \/\ — X 



7(x) = VI -X - - log 7^=, (47) 

Since 7(x) is positive real for x < 1 , in this region 1-instanton amplitude is strongly 
suppresed with order -^e"^ . 
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For the case of ^4 > Tr^m^, both saddle points must be selected because t± steepst 
descent lines are consisitent with original contour and satisfy 



Re[$^(t+)] = Re[$^(t_)]. 



(48) 



Therefore, we get 



dt 

2wt 



;-) 



N 



ie'4 



27rN 



167r-«m4 
A2 



v5 



3_ I 2Tr^m^N .J A \ 

A ^(7^) 



(49) 

Since 7(0;) are pure imaginary in x > 1, the above equation has order A^~^/^ in the 
large N limit. 

Prom the above observation the following truncation about the instanton charge is 
justified. If A satisfies the relation /c^tt^ < A < (k + l)^7r^ {k : positive integer ), we 
truncate the instanton charge as —k <m<k. That is , According to A , the instanton 
charges which are included in the large N limit are changed stepwise as figurel. 

fig.1 



0000 

-o — e — e — o- 



A/pi«2 

Figurel: The instanton charges which must be included in the large N limit (indicated 
by circle ) are changed stepwise according to A. 

This truncation is basically justified in the multi- instanton amplitude because 
become factorized as (cti)^ in the interaction free approximation. But we will see in 
the next section a peculiar dynamics in the multi-instanton interaction make further 
truncation justified in the large N limit. 

Since ai has order N~^/'^ in A > 7r^,we may replace eq (43) to exp(Q;i).But this 
gives order N'^/"^ result for the free energy F = ^ \og[Z/Z^eak]- To see exphcitly that 
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the dilute gas approximation is insufficient we sliow in figure 2 tlie N = 5 result of 
the free energy with the dilute gas approximation and the large N Douglas-Kazakov 
solution which is described later . 

fig.2 



:N=5 
: D-K 



A/piA2 

Figure2: The free energy of A'^ = 5 in the dilute gas approximation and the large N D-K 
solution. 

Next we consider the dilute gas approximation of the Wilson loop. This lead to 



[l7o(Ai, A,) + j^{P,{Ai, ^2)+7i(^i, ^2)) 



sr^N-l iaijA))'' 
2^k=l 



k\ 



iaijA))'' 



n A\A2 

e 2^J^ (50) 



for the Wilson loop. Thus we consider /5i and 71. In the same way as ai ,the large N 
saddle point method leads to /3i,7i = 0{^e^^) for A < n'^ and /9i,7i = 0{^) for 
A > n'^. Here we use the rescaling t2 Nt2 for 71. Thus in the large N hmit the dilute 

gas approximation lead to the value in the weak coupling phase. To sec explicitly the 
insufficiency of the dilute gas approximation we show in figure 3 the N = 5 results of 
the Wilson loop average 



Wi{Ai, A2) = Wi{Ai, A2) - —70(^1,^2)6 



(51) 



in the approximation and compare the results with the large N Boulatov,Daul-Kazakov 
(BDK) solution 



W^''''{A„A2) = W,^^^{A„A2) 



BDK/ 



A,A 



A 



(52) 
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fig.3 (i) fig.3 (ii) 



• : BDK 

+ : N=S 



O 



1 2 3 4 "0 0.2 0.4 0.6 0.8 1 1.2 1.4 

A/pi'^2 A1/pi'^2 

Figures : The Wilson loop average of = 5 in the dilute gas approximation and the 
large BDK solution with (i) Ax = Ajl (ii) A = I.Stt^ 



5 Large N neutrality 

In the preceeding sections we gave the contour integral representation for the multi- 
instanton sector of the partition function and the Wilson loop. The representation has 
a simple gauge group parameter (A^) dependence. That is , 7 have the form such 
as : 

^n^^e-^^W(7W- (53) 

where gii) depends not on for at and weakly depends on for Pk,1k- Then ,If we 
consider the large N limit, the integrals for a, /?, 7 are dominated by the saddle points 
given by the solution of df{t) = 0. Fortunately, it is no need for solving the complicated 
saddle point equations corresponded to the multi-dimensional contour integral in our 
case. The saddle point equations are decoupled for each variables and the solutions of 
the equations are sets of each one dimensional solution which was found by Gross and 
Matytsin 0. Then if we neglect the interaction term M12M23 ■ ■ ■ Mji ( for example aj 
's case ) , aj becomes factorized such as {aiY and the behavior change from exponential 
damping to oscilating at A = vr^m^ observed in the preceeding section still hold for the 
multi-instanton amplitude. But there is a remarkable feature in the interaction term 
of the multi-instanton sector. 

We point out that in the region of A > tt^ special charge configurations of the 
multi-instanton exist. For this purpose let us first consider the following toy contour 
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integral: 



dt 



27rr (t-tc)^' ^^^^ 
where tc is the saddle point of (f){t). We can approximate the above equation in the 
large N limit to a principal value integral along the steepest descent line and a contour 
integral along small semi circle around the saddle point. As the result we obtain finite 

71—1 

value with order N~2~ exp{—N(p{tc)). We call this type of problem "singular saddle 
point problem ". Apparentry, the degree of singularity in ^„ at the saddle point 
determines the power of N. 

In our case such problem occures in a complicated manner. Let's consider the 
singularity of ctfe. The maximal singularity occures when t's are the solution of M^^^ ~ 
0. That is 



-1 




m2 
mi 







m2 







mi 



I 



tl 

t2 



I 



mi 
m2 



-lyv^^ y v-^ / 

The consistency of the above equation requires mi + ■ ■ ■ + = i.e. neutral 
uration. Under this condition the equation has a one parameter solution given 
( t^ \ ( rB±(i \\ mh^ilz^II}! \ 

' ' mi 



^1 

tk-1 
tk 



mi 
rrik 



m2 



{t+l) + 



mi 

mfc_i-| hma 

m2 



irik- 



V 



t 



(55) 
config- 

by 

(56) 



For the case of even k if we consider the following configurations 
(mi, m2, ms, m^, . . . , nik) — (m, — m, m, — m, . . . , — m) 
which we call "neutral" , the above solution simply becomes 

teven t ) ^odd i.^ ~^ -^)' 

On the other hand , under the above charge configuration , the large N saddle 
for any ti are same and given by : 



t 



± 



t+ + t^ + 1^0. 



(57) 

(58) 
points 

(59) 
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When A > rri^n'^, both saddle points must be selected for each variable. Hence the 
saddle points include the following combinations: 

teven = ^± , ^odd = ^ ~(^± + "'-)• (^0) 

In this case the above saddle points ride on the line of maximal singularity and thus 
we encounter the singular saddle point problem. We can also show the "neutral" 
configuration is the only case in which the saddle point completely rides on the line of 
the maximal singularity. Hence the neutral configuration are dominant in the large N 
limit of ak ■ Compared to the one dimensional toy model , we have not the point like 
singularity but the line like singurality and the steepest descent line is also singular 
in this case. This makes a trouble in estimating the power of in a^. But we show 
in the next section numerical calculation indicate that for even k ak has linear scaling 
with respect to N. 

For the case of odd k , the saddle points can not completely ride on the line of 
maximal singularity. But if we weaken the singularity, there is a set of configurations: 

(mi, 1712, ms, 1714,..., rrik-i, ruk) 
— {{m, m, —m, m, ... ,171, —m), (m, —m, —m, m, . . . ,m, —ni), 

(m, —m, m,m, . . . ,m, —m), . . . , (m, — m, m, —m, . . . ,m, m)} (61) 

which we call "next neutral". Under these configurations Mj~^^ = {i = 1, .., j , ..,k) 
are satisfied by the following combinations of saddle points in the region of ^4 > m^vr^ 

{h,...,t^,...,%)^{t^,t^,...,t^). (62) 

Apparently the number of the "next neutral" configurations is k for fixed m. These 
configurations are dominant in the large N limit of 0;^. The numerical calculation in 
the next section showes that a^Ad seems to oscilate with N~^/'^ scahng . 

Next we consider Pk,lk i-e. the connected A;-instanton contribution to the Wilson 
loop. Because we now consider large N limit and the fact that the additive imaginary 
charge in our representation of the Wilson loop has order 1/N , the singular saddle 
point problem occure in the same manner as a^. Thus we conclude that the dominant 
instanton configurations are "neutral" for /3even,7even and "next neutral" for /3odd)7odd- 
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The numerical calculation in the next section show that /3even,7even have also linear 
scaling (N) and (3odd, 7odd have the oscilated iV~^/^ scaling. 

Prom the above argument , we find a peculiar phenomena about the multi-instanton 
configuration in the large N limit. We call the phenomena "large neutrality" . By 
"neutrality" , we would not intend general neutral configuration which satisfy Y,^i i^i — 
0, but intend that all instantons have equal absolute charge and number of positive 
charge instantons and number of negative charge instantons are same or at most differs 
by one. 



6 Numerical calculation 

In this section we show the results of numerical calculation for the connected multi- 
instanton amplitudes ( a's, /3's, 7's) , the free energy and the Wilson loop for finite A^. 
We assume < A < 47r^. Prom the above argument, in this region we can restrict the 

instanton charge as rrii = ±1 in the multi-instanton sector. In the numerical calculation 
we dont perform the integral numerically because wc must deal with highly fluctuated 
integral like the Fourier transformation. But we perform the residue summation which 
is flnite series for flnite in our case. 

First, we consider ak ■ Figl shows the numerical results of ai{A){i = 1, ... ,4) with 
various N and A. We used the following approximation . For acvcn wc restrict the 
"neutral" configuration eq.(56) with m = ±1 and for Oodd restrict the "next neutral" 
configurations eq.(60) with m — ±1 which are dominant in the large limit. 



fig.4 (i) fig.4 (ii) 




Figure4: The connected multi-instanton amplitude aj up to j = A with (i) A'" = 10 and 
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various A/vr^ and (ii) A/vr^ = 2 and various A^. There are remarkable difference between 

eleven 

and Oodd • 



We observe that compared to the a 's in A > vr^ , the a 's in A < tt^ are neghgible. 
Moreover we see remarkable difference between aeven and Oodd with respect to order 
and shape in the region of A > vr^. The shapes of 0:2, ^4 are similar and both have the 
properties of non-oscilating and linear scaling with respect to A^. On the other hand 
ai, as have the properties of oscilating and N'^/"^ scaling . We already observed that ai 
has A^~^/^ scaling and 02 has linear (A^) scaling using the large saddle point method 
0. Hence from the argument in the previous section and these numerical results , we 
conjecture that Oeven have linear (A^) scaling but Oodd have N~^^'^ scaling in the large 
A^ limit. So there is "democlacy" in the interaction between the instantons. That is, 
a2 which is the 2-body force , which is the 4-body force ,. . . are same order with 
respect to A^. 

Figures shows the free energy by our "neutral" multi-instanton scheme 



Z{A) I r dz I 1 



^weak 



[N/2] 



n=l 



(63) 

for A^=3,4 and 5 and the large A^ exact solution of the free energy obtained by Douglas- 
Kazakov. Here we neglect Oodd from the above observation. 
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A = 1.27r2 


A = 1.47r2 


A = I.Qtt^ 


A = l.Sir^ 


A = 27r^ 
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-0.03213 
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-0.08588 



FigureS: The free energy for = 3, 4, 5 approximated by the "neutral" configurations with 
m = ±1 and the large exact solution by Douglas-Kazakov. The neutral multi-instanton 
picture well describe the behavior of the large N D-K solution in the neighborhood of the 
phase transition point. 



The D-K solution is given by the following implicit equations: 



F'{A) 
A 



a 1 

~T ^ 12 • 

b b 

a a 

AK(-)/A 
a 



a 

a 



, , H \ a 1 -) A, 

2' 24 96 ^ 



6^ b 



(64) 
(65) 
(66) 



where F{A) = log Z{A) with the boundary condition -^(vr^) = log Z^^aki'^'^) and 
K{x), E{x) are the complete elliptic integral of first, second kind respectively . The 
above equations are obtained BIPZ-like analysis [n| of the heat kernel representation 
of the partition function. 

The partiton function for A^=3 in our approximation turns into unphysical region 
in A^3.2tt'^ and thus abort in this region. We observe that the graphs for finite 
stand up at y4 ^ vr^ and in the neighborhood of A = vr^ the results for finite seem to 
converge on the large A^ D-K solution as A^ becomes large. But in the region A^2ti'^ 
discrepancy of the finite A^ graphs and the D-K solution becomes large. Since we can 
see that the further approximation with respect to absolute value of instanton charge 
does not make the discrepancy small , we conjecture that the instanton summation 
form of the partition function is a kind of asymptotic series in the region of A > vr^. 

Next we consider the jSkj'jk- We follow the same approximation scheme as the ak- 
That is , for /^even^Teven wc rcstrict the "neutral" configurations with m = ±1 and for 
/^odd^Todd restrict the "next neutral" configurations with m = ±1. 



20 



fig.6 (i-b) 



fig.6 (i-c) 
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fig.6 (iii-b) 
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Figure6: The contribution of the connected multi-instanton ampHtude /3j, ^jto the Wil- 
son loop with (i) iV = 5 , various A = 2Ai (ii)iV = 5,A = 2tP- various A\ (\\\)A = 
27r^, Ai = TT^ , various N are shown 
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FigureG show = 1, . . . , 4) and 7j(j = 0, . . . , 4) with various A^, A, Ai. In fig. 6 we 
observe that compared to the in A > vr^ , the Pj, 7^ in A < tc^ are also neghgible. 
But compared to aj there is no remarkable difi^erence with respect to order between 
/Seven and Podd and between 7even and 7odd in figure6(i)(ii). We should note these are 

= 5 results. And we show in fig.6(iii) both /Seven, 7even have linear scaling with respect 
to A^ contrasted to oscilated damping behavior of Podd, 7odd • Thus in the large A^ limit 
only /3even5 7even are survived and the approximation in which /^odd^Todd are neglected 
fot finite A^ is more nearer to the large A^ exact solution than the approximation in 
which all /3, 7 are included. 

Figure? show the simple Wilson loop average (i.e. winding n = 1 ) for A^ = 3, 4, 5 
by our "neutral" multi-instanton scheme : 



Wi{AuA2) 



O ^ AT M * 



1 X f^-^ /"^'-l 

2^1=0 J 2^{2^k=o 



(67) 

i-i f^^^^A^ + exp(- e2i l^'^'a.^iA))^ 



§ ^P^TT^ exp(- hz^-aUA)) 
where /3odd,7odd are neglected and the large A^ exact solution obtained by Boulatov 
and Daul-Kazakov (BDK)[|Tg]. 
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fig.7 (ii) 
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Figure!: The Wilson loop average with (i) A/" = 3, 4, 5 , various A = 2Ai {n)N = 3, 4, 5, A = 
I.Stt^ various Ai by our "neutral" multi-instanton scheme and for each case showed the 
large N Boulatov,Daul-Kazakov solution. The neutral multi-instanton scheme well de- 
scribes the behavior of the large N BDK solution in the neighborhood of the phase 
transition point 

The large N BDK solution is also given by the following implicit equation : 

W,{A,,A2)= I -^e^i^-^(^). (68) 
Jc Zm 

where f{z) is the resolvent in this system and the contour C encircles the cut oi f{z) 
countercloskwise. Prom the above equation the difference of the phase in the Wilson 
loop average is concentrated on the resolvent . The resolvent is given by the famous 
Wigner's semi-circle law in the weak coupling phase {A < n^): 

A 7 A I 

m ^ ^ - -^^/z'-A/A (69) 

On the other hand in the strong coupling phase the resolvent is given by 

A 7 2 I h 

f{z) ^ — + - J{a^-z^W-z^)7r{--, -) (70) 

^ QjZ Z Qj 

where a, h are the same parameter as the equations of the free energy and 7r(c, x) is the 
complete elliptic integral of third kind. 
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We observe that the graphs for finite stand up at A ~ vr^ and in the neighborhood 
oi A = n"^ the resuhs for finite N seem to converge on the large BDK solution as 
becomes large. But we also observe that the disagreement becomes large in the region 
A > 27r^ which is same as the free energy. Hence the instanton summation form of the 
Wilson loop is also a kind of asymptotic series and divergent in the region of A > vr^. 

Hence we conclude that viewed from weak coupling phase the large N phase tran- 
sition is physically caused by the "neutral" multi-instanton. 



7 conclusion 

We have analyzed the multi-instanton amplitudes in terms of the contour integral 
representation and found that the represenation make clear the properties and the role 
of the multi-instanton in the large phase transition. In particular the "neutral" 
configurations of the even number instantons are essential. We also showed that our 
"neutral" multi-instanton scheme well describes the behavior of the free energy and 
the Wilson loop around the critical point Ac = vr^ up to A ~ 27r^. And we show that 
is imposibble in the dilute gas approximation. 

But we observe the agreement region of our results and the large A^ exact solution 
tends to become narrow as A^ becomes large. In addition we can see that the further 
approximation with respect to absolute value of instanton charge does not make the 
discrepancy small. So we conjecture the form of the instanton summation is a kind of 
asymptotic series. For example , assume that S{z) is defined nonperturbatively and has 
an asymptotic expansion S{z) =^ Y^'k'=o^kz''- By truncating the series with apropriate 
region of summation in the neighborhood of 2; = , we can approach to the true value. 
The same scenario works for our case of the "neutral" multi-instanton truncation. 
Thus we conclude that our "neutral" multi-instanton scheme is a good truncation in 
the strong coupling region up to A ~ 27r^ at large A^. 



Recently the field strength correlators in the large A^ limit are determined [|l^] using 



the abelianization technique for the path integral |T^. To explore whether our method 



can describe the behavior of these correlators is our next task. 
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